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ABSTRACT 

Aims. The intracluster light (ICL) is a faint diffuse stellar component in clusters made of stars not bound to individual galaxies. We 
have carried out a large scale study of this component in the nearby Virgo cluster. 

Methods. The diffuse light is traced using planetary nebulae (PNe). The surveyed areas were observed with a narrow-band filter 
centered on the [OIII]/i5007 A emission line redshifted to the Virgo cluster distance (the on-band image), and a broad-band filter (the 
off-band image). For some fields, additional narrow band imaging data corresponding to the Ha emission were also obtained. The 
PNe are detected in the on-band image due to their strong emission in the [OIII]/i5007 A line, but disappear in the off-band image. 
The contribution of Ly-o- emitters at z = 3.14 are corrected statistically using blank field surveys, when the Ha image at the field 
position is not available. 

Results. We have surveyed a total area of 3.3 square degrees in the Virgo cluster with eleven fields located at different radial distances. 
Those fields located at smaller radii than 80 arcmin from the cluster center contain most of the detected diffuse light. In this central 
region of the cluster, the ICL has a surface brightness in the range p. B - 28.8 - 30 mag arsec~ 2 , it is not uniformly distributed, and 
represents about 7% of the total galaxy light in this area. At distances larger than 80 arcmin the ICL is confined to single fields and 
individual sub-structures, e.g. in the sub-clump B, the M60/M59 group. For several fields at 2 and 3 degrees from the Virgo cluster 
center we set only upper limits. 

Conclusions. These results indicate that the ICL is not homogeneously distributed in the Virgo core, and it is concentrated in the high 
density regions of the Virgo cluster, e.g. the cluster core and other sub-structures. Outside these regions, the ICL is confined within 
areas of ~ 100 kpc in size, where tidal effects may be at work. These observational results link the formation of the ICL with the 
formation history of the most luminous cluster galaxies. 

Key words, galaxies: clusters: individual (Virgo cluster) - galaxies: halos - planetary nebulae: general 



■ 1. Introduction importan t transformatio ns over larger timescales due to galaxy 

)\ ' mergers (Mihoi, I2004 and references therein). According to 

H ; Galaxy clusters are the most massive gravitationaUy bound these resultS; , evolution is not frozen at the earliest 

structures known in the Universe. In a volume of a few Mpc 3 Qf ±& fOTmation> afj described in earlier works 

they contain a large number of galaxies, which interact with each (Richstone & Malumuthi 1983; MerritiCHi, but is still active 

other, the hot gas and with the global cluster gravitational poten- tof i av 

tial. Thus, in these high density regions the evolution of galaxies 

may be influenced by the environment. ^ Hierarchical structure formation theories predict that mas- 

With new numerical algorithms and fast computers, we Slv , e § alax y clusters are built through the infall of matter, i.e. 

are able to simulate large volumes of the Universe and fol- ? alaxies ,; g°T rr^ T l'brS ^ 6 t§ 

low the evolution of galaxies in different environments as jWhtte & Reesl pZl | West et al l [lpl). Since mass accretion 

function of redshift. These simulations show that galaxies un- 18 stl " ac f u e toda y< nearb y 8 alax y J clu u sters m ^ b „ e at dlfferent 

dergo important transformations over time in cluster environ- e P och J s of their u mas s assembl y< and they may still be in a un- 

ments. They can evolve over short timescales (less than 1 mlxed sta f The dynamical state of a galaxy cluster has been 

Gyr) due to different mechanisms such as fast interactions traditionally measured by determining the degree of substruc- 

between galaxies and with the cluster gravitation al potentia l ture f ln the 8 alax y distribution and in the X-ray emitting hot in- 



(Merritt, 1984; Moore et al., 1996), gas stripping (Gunn & Gott, traclu f er medium, but the inferred results do not always agree. 
1973 lOuiks et all 120001) and starvation (lLarson et all 11980: (see| Ramella et al| | 2007i and references therein) 



Bekki et all 12002b . Moreover, in the cluster central region, The gravitational forces acting on galaxies during the cluster 



where the galaxy density is high, galaxies can also undergo formation and its evolution unbind a fraction of their stars, which 

then end-up orbiting in the intracluster region. This unbound dif- 

Send offprint requests to: N. Castro-Rodriguez fuse stellar component is the so-called intracluster light (ICL). 
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Table 1. Summary of the field positions, filter characteristics and exposure times for the narrow-band (on-band) images. 



[OIII] Ha 



Field 


RA(J2000) 


DEC (J2000) 


A c 


FWHM 


Exposure 


Sfwhm 


A c 


FWHM 


Exposure 


Sfwhm 




(hh mm ss) 


(° ' ") 


(A) 


(A) 


(s) 


(arcsec) 


(A) 


(A) 


(s) 


(arcsec) 


Core 


12 27 48 


+ 13 18 46 


5023 


80 


24000 


1.20 


— 


— 


— 


— 


FCJ 


12 30 39 


+ 12 38 10 


5027 


44 


18000 


1.40 


— 


— 


— 


— 


LPC 


12 25 32 


+ 12 14 39 


5027 


60 


27000 


1.40 










SUBC 


12 25 47 


+ 12 43 58 


5021 


74 


3600 


0.68 


6607 


101 


8728 


0.70 


LPRX 


12 28 00 


+09 30 00 


5027 


60 


27000 


1.80 










LPE 


12 17 27 


+ 13 41 30 


5027 


60 


27000 


1.30 










LPS 


12 27 34 


+10 43 42 


5027 


60 


27000 


1.30 










LSF 


12 38 33 


+ 12 08 41 


5023 


80 


33000 


1.08 










RCN1 


12 26 13 


+ 14 08 03 


5023 


80 


24000 


1.20 










SUB2 


12 24 13 


+ 14 09 01 


5021 


74 


3600 


1.00 


6607 


101 


9000 


1.10 


SUB3 


12 24 32 


+ 15 21 21 


5021 


74 


3600 


1.00 


6607 


101 


9000 


1.05 



The fields listed above the horizontal line are located within 80 arcmin of the Virgo cluster center; the data reduction is described in previous 
work (FCJ: Arnaboldi et al. 2002; SUBC: Arnaboldi et al. 2003; Core and LPC: Aguerri et al. 2005). The fields listed below the horizontal 
line are at distances between 80 and 215 arcmin of Virgo cluster center. 



Table 2. Summary of the filter characteristics and exposure times for the broad-band (off-band) images in the various fields. 









B-band 








V-band 








R-band 




Field 


Ac 


FWHM 


Exposure 


Sfwhm 


Ac 


FWHM 


Exposure 


Sfwhm 


Ac 


FWHM 


Exposure 


Sfwhm 




(A) 


(A) 


(s) 


(arcsec) 


(A) 


(A) 


(s) 


(arcsec) 


(A) 


(A) 


(s) 


(arcsec) 


Core 










5395 


894 


2400 


1.20 










FCJ 










5300 


267 


6000 


1.4 










LPC 


4407 


1022 


5400 


1.50 


















SUBC 










5500 


956 


900 


1.00 


6500 


1130 


720 


0.86 


LPRX 


4407 


1022 


15600 


1.60 


















LPE 


4407 


1022 


15600 


1.50 


















LPS 


4407 


1022 


15600 


1.50 


















LSF 










5395 


894 


14300 


1.10 










RCN1 










5395 


894 


2400 


1.20 










SUB2 










5500 


956 


3600 


1.10 


6500 


1130 


1440 


0.80 


SUB3 










5500 


956 


3600 


0.80 


6500 


1130 


1440 


0.90 



Numerical simulations show that the amount of ICL in clusters 
depends on the cluste r mass a nd dynamical state {Murante et al.L 
12004 IWillman et all 12004 ISommer-Larsen et all 120051) . A 
more massive, older, dynamically evolved cluster may con- 
tain a larger amount of ICL than a less massive or a dynami - 
cally younger system dRudick et all 120061: IMurante et alll2007l) . 
Furthermore, in a highly evolved dynamically old cluster the ICL 
morphology would be more diffused, with relative few streams, 
while a dynamically young cluster dominated by groups still in 
the proce ss of merging is like ly to be dominated by ICL found in 
streams (Rudi cket al.L 12009). Therefore the study of the amount, 
distribution and kinematics of the ICL may provide information 
on the cluster accretion history and evolutionary state, as well as 
about the evolution of cluster galaxies. 



The ICL 
techn iq ues, e . 
2002, 12004" 



has been mapped using several observational 
deep broad ba nd im aging 
iMinosetall 12001 
Krick & Bernst ein. 2007), and the detection of individual 



r 



_ (Feldmeieret al 
Zibe tti et all 12005 



stars associated wi th the diffuse s t ellar c omponent, both red 
giant branch stars jFerguson et all 119981; iDurrell et all [2002; 



Williams et al., 



Arnaboldi et al 



2007]) and intracluster planetary nebulae (ICPNe 
19961120021 12003b iFeldmeier et al.L[l998ll2003l 
2004at lAguerri et all 120051 l2006h . In particular, ICPNe are the 



only componen t of the ICL whose kinematics c a n be measured 
dFreeman et all l2000t lArnaboldi et all 120041: iGerhard et ail 
l2005Ll2007tlDohertv et alll2009l) . 



Several studies investigate the properties of the ICL in the 
core of the nearby Virgo cluster. E xpanding on our earlier imag- 
ing work in the Virgo cluster core (Arnaboldi et all 120021 [2003; 
Oka mura et all 120021 ; lAguerri et aL , 2005b . we have now com- 
pleted a survey campaign of the ICL distribution on larger scales, 
outside the center of the Virgo cluster. In total, we have now 
covered more than 3 square degrees in Virgo, at eleven different 
positions in the cluster and at distances between 80 arcmin and 
some 100 arcmin from the Virgo cluster center. In several of the 
new fields, the ICL is at least two magnitudes fainter than in the 
core region. These new results are in agreement w i th obs erva- 
tions of intermediate redshift clusters (Zibetti et al., 2005), and 
with the results of recent high resolution hydrod ynamical simu- 
lations of cluster formation in a A CDM universe (IMurante et all 
12004 iRudick et all 120061; IMurante et ail 120071) . which predict 
that the ICL is more centrally concentrated than cluster galax- 
ies and that the largest portion of the ICL is formed during the 
assembly of the most luminous cluster galaxies. 

In the following, we will present the results of a systematic 
and homogeneous ICPNe survey in regions outside the Virgo 
core and discuss the results jointly with ICPN survey data avail- 
able for the Virgo cluster. In Section 2 we will present the data 
for the new fields in the Virgo cluster. The technique for the se- 
lection of ICPNe will be summarized in Section 3. The possi- 
ble contamination of the ICPN catalog by background emission- 
line objects is quantified in Section 4. In Section 5 we use the 
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number of detected ICPNe to infer the surface brightness of the 
light at the fields' position and compare our results with pre- 
vious works. We discuss the distribution of the ICL from the 
Virgo core to the outer distant fields at ~ 200' = 0.8 Mpc, and 
its implications for the current models for the formation of ICL 
in clusters, in Section 6. In Section 7 we set limits on the con- 
tribution from a possible homogeneous luminous density in the 
nearby Universe to our PN sample, and finally summarize our 
conclusions in Section 8. 

In what follows we con sider the center of th e Virgo cluster 
at the position of Cluster A Bingg elTet alj dl987l) . at the peak of 
the luminosity density of all galaxies in the Virgo cluster, e.g. 
a(72000) = 12 : 27 : 49.95 5(72000) = 13 : 01 : 24.57. 
We adopt a distance of 15 Mpc to the Virgo cluster center, thus 
l"=72.7pc. 

2. Observations and Data Reduction 

Previous ICPN surveys i n the Virgo cluster were car r ied out 
mostly in the V i rgo core (iFeldmeier et al.L [19981 12003L l2004d: 
lArnaboldi etaiT. l2002l l2003b lAguerri et all 120051) and around 
M49 (IFeldmeier et alll2*004al) . 

We have now collected additional data in several fields lo- 
cated at larger radial distances from the center of the Virgo clus- 
ter. In total we h ave now six new fie lds, which in addition to 
those surveyed in Agu erri et al.l d2005l) cover a total area of 3.3 
square degrees. The Virgo cluster SDSS image and the location 
of the eleven fields (five fields surveyed in previous studies by 
our group and the new six ones) are shown in Figure 1 . We in- 
dicate also the areas surveyed by Feldmei er et al.l (1 19981 2003, 
2004a) and the name of the brightest Virgo galaxies. The results 
of the survey for ICPNe in the LPRX, LPE, LPS, LSF, SUB2, 
SUB3 field are presented here for the first time. In Table 1 and 2 
we give the eleven fields' IDs and locations. 

Three of the new fields were observed during two observ- 
ing runs in March 2002 and 2004 with the Wide Field Camera 
(WFC) on the 2.5 m Isaac Newton Telescope (INT) at Roque de 
los Muchachos Observatory (La Palma, see Tables 1 and 2). We 
refer to them hereafter as the LPRX, LPE and LPS fields. The 
images from this camera consists of a mosaic of four CCDs with 
a total field of view of 34' x 34'. The pixel scale is 0". 3 3 3, and 
the mean readout noise and gain of the detectors are 6.1 ADU 
pixel 1 and 2.8 e~ ADU 1 , respectively. We imaged these fields 
through a 60A wide filter centred at 5027 A, which contains the 
wavelength of the [OIII]/15007A emission at the redshift of the 
Virgo cluster. We will refer to this narrow-band filter as the "on- 
band" filter (see Table 1). We also took an additional broad-band 
image per field through a filter centred on the B band (4407 A) 
and 1022 A wide. The broad-band filter is referred to as the "off- 
band" (see Table 2). The total exposures were 27000 sec for the 
on-band images and 15600 sec for the off-band. 

In April 2004 we acquired another field (LSF) with the Wide 
Field Imager (WFI) on the ESO/MPI 2.2 m telescope located at 
La Silla Observatory (Chile). The full image from this camera 
consists of a mosaic of eight CCD images, covering 34' x 34' on 
the sky. Each CCD has a pixel scale of 0".238, and an average 
readout noise and gain of 4.5 ADU pixel 1 and 2.2 e~ ADU 1 , 
respectively. The field was imaged through an on-band 80A wide 
filter, centred at 5023A. In addition we also imaged this field 
through an off-band V-band 894A wide filter, centred at 5395A 
(see Table 1 and 2). 

From January to April 2004 two fields in the Virgo Cluster 
(SUB2 and SUB3) were imaged with the Suprime-Cam 10kx8k 



mosaic camera, at the p rime focus of the 8.2 m Subaru telescope 
(MiyazakieFm|2002|). The field of view of these images covers 
an area of 34'x27', with a pixel size of 0".2. The fields were im- 
aged through two on-band filters and two off-band, the standard 
V and R filters (see Tables 1 and 2). 

The new fields at larger cluster distances from M87 are com- 
plementary to those studied in the Virgo cluster core; their pa- 
rameters are reported in Table 1 and 2. The data reduction and 
calibration of the new fields, LPRX, LPE, LPS and LSF, fol- 
low the procedure described in Arnabo ldi et al.l (2002, hereafter 
Paper I). The data reduction is performed using the MSCRED 
package in IRAfQ- Calibrations are done using several Landolt 
fields for the off-band filters, and spectrophotometric stars for 
the on-bands. Fluxes are normalized to the AB magnitude sys- 
tem, following [fh euns & Warrenl d 1997b . A detailed description 
of the calibrations steps and the relation between the A B mag- 
nitude s and the m(5007) [OIII] magnitude introduced bv lJacobvl 
(1989) is given in Paper I. Table 3 shows the conversion constant 
(K) between the two magnitude systems: m(5007)=myiB+K. 

The data for the additional new fields observed with the 
SuprimeCam at the SUBARU telescope (SUB2 and SUB3) 
are processed using the standard data reduction package devel- 
oped for SuprimeCam data. The photometric calibrations for the 
SUB2 and SUB3 data are obtained using Landolt fields for the 
broad bands, and spectrophotometric stars for the narrow band0 

For completeness, Table 1 and 2 provide a summary of the 
fields positions, filter characteristics and exposure times for on- 
band and off band exposures for the surveyed areas in the Virgo 
core. Throughout this work we will use the notation m„ and nib 
to refer the on-band and off-band magnitudes of the objects in 
the AB photometric system. 

3. Detection and selection of ICPN candidates 

In this work, the diffuse light in a Virgo cluster field is mea- 
sured by the number density of detected PNe. The PNe are iden- 
tified via their strong emission in the [OIII]/15007A line and, 
in the two SUB2, SUB3 fields, may be confirmed by the ad- 
ditional detection of the associated Ha emission. The obser- 
vational techniques used for the PNe detection have been suc- 
cessfully applied in the Virg o cluster, in nearby groups (Leo; 
Castro-Rodriguez et al., 2003), and in compact groups (HCG44; 
Agu erri et all 120061) . In a few Virgo fields we carried out 
the spectroscopic follow-up, resulting in a large fraction (up 
to 80% for some fields) of spectroscopica lly confirmed PNe 
(lArnaboldi et all 12004 iDohertv et all 120091) . 



3.1. [OIII] "on-off" band survey of ICPNe - Catalog selection 
with the color magnitude diagram (CMD) 

The ICPN candidates in the fields LPRX, LPE, LPS and LSF 
were identified using the on-off band technique, pioneered by 
Ciardullo et al . (198 9) and automatized for large areas by 
Arnabold i et al.1 (120021 Pa per I), using selectio n criteria based on 
color excess discussed by Theuns & Warrery d 1997b . On the ba- 
sis of the strong line emission of a PN at the [OIII]/15007A line, 
the photometric ICPN candidates can be identified as point-like 
sources detected in the [OIII] image, the on-band and, because 
of their faint continuum, they cannot be detected in the off-band 



1 For a detailed disc ussion of the mosaic data reduction, we refer the 
reader to lAlcala etal] j2002h . 

2 Zero points for the Subaru data are zo.Ha = 24.0, zo.iom] = 24.6, and 
Zo,v+« = 27.8 for a 1 sec exposure 
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Fig. 1. Finding chart of the Virgo fie lds discussed in this p aper. The circles indicate radial distances of 1, 2, 3 and 3.5 degrees from 
the Virgo cluster center identified bv lBinggeli eTail dl987l) : a(72000) = 12 : 27 : 49.95 (5(72000) = 13:01: 24.57. Open squares 
indicate the field positions whose data were analyzed by our team. The area covered by each field is : 16x16 arcmin 2 for the FCJ 
field, 34x27 arcmin 2 for SUBC, SUB1 and SUB2 fields. Th e remaining field s have areas of 34x3 4 arcmin 2 (Core, RCN1, LSF and 
LP* fields). The shaded fields indicate the areas surveyed bv lFeldmeier et al] dl99 8, 2003, 2004a); they are designated as F04-n in 
the text, where n is the number given in the respective shaded area on the figure. 



image. The on-band filter is a few tens of A wide filter (60 or 80 
A, see Table 1 and 2) centred at the wavelength of the redshifted 
[OIII]5007A emission of the PN population at a given distance. 
In Paper I, this technique was standardized and automatized to 
select ICPNe in wide field images. We give a brief summary 
here. 

We detect all objects on the on-band imag es. We then per- 
form aperture photometry with SExtractor dBertin & Arnoutsl 
1996) on the on-band image at the sources' position, and in 
the corresponding aperture in the off-band image. All point-like 
sources detected on the on-band image are placed in a color- 
magnitude diagram (CMD), m„-m,i, vs m„, and are classified ac- 
cording to positions in this diagram. The ICPN candidates are 
point-like objects with no detected continuum emission and ob- 
served EW greater than the color excess curve corresponding to 
EW=100A and convolved with the photometric error, see Paper 
I for more details. Figure|2]shows the CMDs for the LPRX, LPE, 
LPS, and LSF fields; the selected ICPN candidates are indicated 
by asterisk and square symbols. 



Table 3. Limiting magnitudes of the fields in the on-band, off- 
band filter, and the K constant 



Field 


m,„„(5007) 


m/,/,,,, 


K 


Core 


27.21 


24.75 


2.51 


FCJ 


27.01 


24.58 


2.50 


LPC 


27.52 


25.40 


3.02 


SUBC 


28.10 




2.49 


LPRX 


27.49 


25.05 


3.02 


LPE 


27.83 


25.93 


3.02 


LPS 


27.69 


25.91 


3.02 


LSF 


28.20 


26.32 


2.51 


RCN1 


26.71 


25.63 


2.51 


SUB2 


27.99 


26.81 


2.49 


SUB3 


27.74 


26.65 


2.49 



The magnitude limit of our emission-line catalog is evalu- 
ated using Monte Carlo simulations. We distribute a large sam- 
ple of point-like objects (~ 10000) randomly on the scientific 
images. We then perform the photometry of the objects using 
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Fig. 2. Color-magnitude diagrams (CMDs) for all sources in the LPRX, LPE, LPS and LSF fields. The horizontal lines indicate 
the color excess of emission line objects with an observed EW=100A. The inclined continuous lines show the ;«/, magnitude 
corresponding to 1.0 cr above the sky, in the off-band images. The curved lines limit the regions where 99% and 99.9% of the 
simulated continuum objects would fall in this diagram, given the photometric errors. The dashed lines limit those regions where 
84% and 97.5% of the simulated objects with color excess relative to an observed EW=100A would fall in this diagram, as function 
of magnitude and the photometric errors. The asterisks represent point-like sources with off-band magnitudes fainter than the 
detection limit in the off-band image as measured by SExtractor; the squares represent point-like objects with an off-band emission 
fainter than the limiting magnitude in the off-band image. 



the same procedure as for the real emission sources. The limit- 
ing magnitude is then defined as the magnitude at which 50% 
of the simulated population is retrieved. In Table 3 we list the 
limiting magnitudes for the on-band image in the m(5007) pho- 
tometric system, and for the off-band magnitude in the Johnson 
system. We list also the conversion constant K between the AB- 
system and m(5007). According to Table 3 the "shallowest" field 
is RCN1 (m /im (5007) = 26.71) and the deepest field is LSF 
(m, /m (5007) = 28.20), with m /iffl (5007) fainter than 27.0 for all 
the other survey fields. 

At the distance of 15 Mpc, the bright cut-off of the PN lu- 
minosity function (PNLF) falls at the apparent magnitude of 
m(5007)= 26.383 which is brighter than the limiting magnitude 
in all surveyed fields. 

On the basis of the CMD selection, we detect 9, 20, 7 and 
23 emission-line point-like objects brighter than the m;,„,(5007) 
in LPRX, LPE, LPS, and LSF, respectively. These emission 

3 This is computed us ing the absolute magnit ude for the PNLF bright 
cut-off of M; 007 = -4.5 dCiardullo et alll20023) and a distance modulus 
of 30.88. 



sources include ICPN candidates and background galaxies. In 
Appendix A we give the coordinates and m(5007) magnitudes 
of these sources. 

3.2. [Olll] and Ha survey of ICPNe - Catalog selection with 
the 2-color diagram (2-CD) 

A more secure selection of the ICPN candidates is possible in the 
in the SUB2 and SUB3 fields, using the additional information 
from the Ha images which eliminates the contamination from 
background galaxies (Okamura et al ., 200 2:lMihos et all 12009). 
The selection procedure is described in lArnaboldi et alj (2003, 
hereafter Paper II), on the basis of the instrumental magnitudes 
measured for sources in the fields, normalized to 1 sec exposure. 

We use SExtractor for the detection of point-like objects in 
SUB2 and SUB3 [OIII] on band images. Once the objects are 
extracted, we measure fluxes at the sources' positions in the 
Ha, V and R band images in matching apertures. We then use 
a two-color diagram (2-CD), [OIII]-Ha vs. [OIII]-(V+R), for 
the selection of the ICPN candidates. In this 2-CD, continuum 
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sources, single-li ne and two-line emi tters, either PNe or com- 
pact HII regions dGerhard et al.L 12002b . occupy different regions 
of the plane. We establish some ad-hoc selection rules for this 
specific set of images, and the detection/photometry described 
above. Simulated population of PNe, single line emitters, and 
continuum sources are constructed in the 2-CD, and are used 
to outline the regions inhabited by the different kind of objects. 
Differently from Paper II, we have the zero points of the photo- 
metric calibration for the narrow bands, [OIII] and Ha, and for 
the broad band images, V and R, and the selection criteria for the 
different classes of objects translate into these color constraints: 

- The simulated PN population consists of point-like objects 
with a color excess in [OIII], e.g. m n - m^v+R) < -1.3, 
and have line r atios \ OIII]/(Ha + [Nil]) larger than 3 
dCiardullo et all |2002a), because they populate the brighter 
part of the PNLF. This line ratio translates into an [OIII] vs. 
Ha color excess of m„ - m.n a < -1. Therefore the ICPN- 
like candidates are point-like sources detected in both [OIII] 
and Ha images, whose colors are m„ - m^v+R) < -1-3 and 
m n — niHa < -1; the simulated objects are plotted as open 
squares in Fig. [3] . 

- The continuum sources are simulated with a flat spectral en- 
ergy distribution which translates into m n - mu a - ni n - 
m(v+R) = colors; these are plotted as open triangles in 
Fig.E] 

- Single line emitters in [OIII] are those objects detected in 
the [OIII] image, and have upper limits in either Ha or in 
the V+R continuum image, or have no detections in those 
images. These are plotted as asterisks in Fig. [3] Bright m„ 
sources without Ha emission are background, large EW Ly- 
a emitters, while fainter m„ sources without Ha emission can 
be either PNe or background objects. Because of the different 
detection limits and noise in the Ha and V+R image, they are 
scattered in the different regions of the 2-CD. See Paper II for 
additional details. 

- Compact HII regions are sources with faint continuum, m n - 
niHa > and m„ - m^y+R) < 0, because of the stronger Ha 
line emission. 

The different class of simulated objects are added to the sci- 
entific images and then recovered using the same automatic pro- 
cedure as for the real objects. Fig|3]shows the 2-CDs of the sim- 
ulated objects for the SUB2 and SUB3 fields. 

The limiting magnitude of the [OIII], Ha, V and R band im- 
ages are listed in Table 3. The m/,,„(5007) for the on-band [OIII] 
images of the SUB2 and SUB3 fields is 1.5 magnitude fainter 
than the apparent magnitude of the PNLF bright cut-off at 15 
Mpc distance. The 2-CDs and CMDs of the emission objects in 
SUB2 and SUB3 images are shown in Fig. [4] 

3.3. Comparison between the two selection procedures 

For the SUB2/SUB3 fields, we can compare the two methods for 
detecting ICPNe either with CMD and 2-CDs. Fig. [4] shows the 
CMD for these two fields: we have obtained 13 and 11 objects 
in SUB2 and SUB3 brighter than the m/, m (5007) in these fields, 
that are classified as ICPNe on the basis of a color excess criteria 
in the CMD. Fig. [4] also shows the distribution of those objects 
on the 2-CDs. It turns out that most of the objects classified as 
ICPNe in the CMDs are located in the region occupied by sin- 
gle line emitter in the 2-CDs, i.e. they do have a detectable Ha 
emission at the position of the [OIII] source on the Ha image. 
Four objects, one in SUB2 and three in SUB3, which are ICPN 
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Fig. 3. Color-color diagrams (m„ -niHa) vs.( m„ -OT(v + r)) of sim- 
ulated objects for SUB2 (a) and SUB3 (b). The different symbols 
represent: continuum objects (open triangles), single line emit- 
ters (asterisks), and PNe (open squares). Vertical and horizontal 
dashed lines show objects with m n -mu a = m„-m^v+R) — 0. The 
vertical dashed-dotted line represents objects with m n -m^v+R) = 
1.3. The horizontal dotted and dash-dotted lines show objects 
with m„ - rn.na = -1 and -2 for SUB2 and SUB3. The con- 
tinuum and/or Ha magnitudes of the simulated objects with no 
continuum and/or no Ha detection were set arbitrarily to 40. 



candidates according to the CMDs, are also two-line emitters 
according to the 2-CD. 

From the 2-CDs, we conclude that we detect one ICPN in 
SUB2 and three ICPNe in SUB3 brighter than the limiting mag- 
nitudes in [OIII] and Ha; the four ICPNe are indicated by open 
diamonds in the right panel of Fig. [4] Furthermore, no HII com- 
pact regions are detected in these two fields. 

In empty fields, the sources selected on the basis of 
a color excess in a narrow band, e.g. in the 4500 - 
5000 A wavelength interval, vs. a broad band filter se- 
lected catalog, are most probably Ly-a background galax- 
ies at redshift ^ 3.1. The identification of these sources 
as Ly-a emitters has been confirmed with spectroscopic 
follow-up dKudritzki et"aU 120001; ICastro-Rodriguez et aklllOOH; 
iGawiser et al.Ll2(XMlRaucIiet al., 2008). Therefore when only a 



single deep [OIII] image is available to produce the catalog of 
point-like emission line objects for a given Virgo field, we need 
to subtract statistically the background population of the Ly-a 
emitters to select the over-density of objects associated with the 
PN population. 




Fig. 4. CMD (left panels) and 2-CD (right panels) of the sources in the SUB2 field (top panels) and SUB3 field (bottom panels) 
images. On the left plots: the squares and asterisks represent objects classified as ICPNe with CMDs. On the right: all CMD- 
identified ICPNe are indicated with open squares; open diamonds show the ICPNe selected from the 2-CD. The continuum and/or 
Ha magnitudes of those objects with no-continuum and/or Ha detections were assigned arbitrarily, as in Fig. [3] 



We now discuss the statistical subtraction of such back- 
ground population in the Virgo cluster fields in Section |4] 

4. On the nature of the emission-line sources in the 
Virgo outer fields 

In addition to a population of ICPNe at the Virgo cluster dis- 
tance, in the CMD-selected source list with color excess, we 
expect a contribution from misclassified faint continuum ob- 
jects, [Oil] emitters at intermediate redshift and Ly-a back- 
ground galaxies at z« 3.1. We discuss the contribution from faint 
continuum sources in Section 14.11 and background galaxies in 
Section [4~2l with the goal of selecting the most likely PN popu- 
lation at these field positions. 

4.1. Contamination by faint continuum objects 

It is possible that our catalogs contain misclassified faint con- 
tinuum objects. Because of the photometric errors and a steeply 
rising luminosity function, a non-negligible number of faint con- 
tinuum objects, which would have m„ fainter than the limiting 
magnitude in the on-band image, may end up having a measured 
m„ brighter than m/ im (5007). As the same time, due to the photo- 
metric errors, the measured off-band magnitude may be fainter 
than what is to be expected, resulting in a color excess for these 



continuum sources as for the ICPNe. Given the large number of 
stars in the field, the number of such sources can be significant, 
if the off-band is not sufficiently deep. 

We named this as the "spillover" effect in Paper III, and we 
refer to it for additional details. The number of continuum ob- 
jects that "spilled over" can be estimated from simulations of 
point-like continuum objects with a LF extrapolated to 3 mag- 
nitudes fainter than m/,„,(5007) and distributed randomly on the 
scientific frames. We then detect the simulated objects with the 
same criteria as for the real ones, and determine the number 
of simulated objects brighter than ra/ im (5007) in regions of the 
CMD where the ICPN photometric candidates are located. After 
scaling the simulations to the number of observed continuum 
objects, the final numbers of spilled-over contaminants brighter 
than m/„„(5007) of the image are found to be zero in all fields. 
This result comes about because of the very deep off-band im- 
ages acquired for these fields (see the discussion in Paper III). 

4.2. Contribution from background galaxies 

Background emission-line galaxies are detected in [OIIIJ 
on-band photometric s urvey s dCiard ullo et all |2002b; 
ICastro-Rodrfguez et all l2003h iGronwall etaU I2007I) . 
According to the central wavelength and limiting flux of 
the on-band images, the surveys are sensitive to the [Oil] 3727 A 
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Fig. 5. The luminosity function (LF) of the emission-line obje cts detected in LPRX, LPE, LPS, LSF, SUB 2 and SUB3 is indicated 
by large full red squares. For comparison, the Ly-a LF from lKudritzki et ail d2000h (open green circles), Ciardul lo et alj (|2002b) 
(open blue triangles) and ICastro-R odrfgue z et alj d2003l) (the Leo field; fille d orange circles) are a lso reported in each panel. The 
continuum line shows the luminosity function of the Ly-a candidates from iGronwall et all d2007l) . The LFs of the Ly-a emitters 
were scaled to the effective surveyed area and wavelength range of the different observed fields. 



emission line of star-burst galaxies at z=0.35 and the Ly-a line 
of young galaxies at z=3.1. 

[Oil] emission line at z ~ 0.35: We do not expect a sig- 
nificant contamination by [Oil] emitters at medium redshifts in 
the final source catalogs because the emission-line objects have 
a color excess which implies observed EWs greater than 100A. 
Previous studies of [Oil] emitters at z=0.35 have not found ob- 
jects with observed EWs greater than 9 5A dColless et all 1 19901: 
Hammer et aflll997l:lHo"gg et al.Lll998l) . The long-slit search for 
low surf ace brightness Ly-a with fluxes of a fewxl0~ 18 erg s _1 
cm -2 bv lRauch et al.l d2008l) also reports a small number density 
of [Oil] emitters at these faint fluxes. 

Ly-a line emitters at z-3.1: On-band/off-band photometric 
surveys of bl ank fields provide a non-negligible density of emis- 
sion objects dCiardullo et aUl2002bt IGronwall et al.Ll2007l) . The 
spectroscopic follow-up of these e mission line objects s howed 
that they are Ly-a- emitters at z=3.1 dGawiser et al I l2007h . 

In what follows we adopt a statistical approach to cor- 
rect for the contribution of the z=3.1 Ly-a background ob- 
jects, by comparing the number density of the "PN-like" 
objects in the outer Virgo fields with those measured in 
"blank-fields" dKudritzki et all 12 000; Cia rdullo et all l2002bt 
ICastro-Rodrfguez et all 120031: IGronwall et all 120071) . normal- 
ized by the effective sampled volume. 



4.3. Are ICPNe detected in the Virgo fields outside the 
cluster core? 



We thus estimate the fraction of background galaxies present 
in the current catalogs of point-like sources by comparing 
their LF and number densities with those of Ly-a emitters de- 
tected in similar [OIII] on band surveys for objects brighter 
than m(5007)= 27.0 and normalized to the effective sam- 
pled volumes. We shall compare the different sample for mag- 
nitudes brighter than m(5007)= 27.0; at fainter magnitudes, 
the automatic classification of pont-like(unresolved) vs. ex- 
tended(resolved) based on total magnitude vs. small central aper- 
ture magnitude plot (see Paper I) becomes unreliable and not 
uniform among the fields' ICPN catalogs and the comparison 
blank field catalogs. 

Fig. [5] compares the LFs of the point-like emission-line ob- 
jects detected in the LPRX, LPE, LPS, LSF, SUB2 and SUB3 
fields with the LF of the Ly-a population from different sur- 
veys. We then evaluate the number of expected Ly-a emitters in 
our fie l d brighter than m(5007)= 27.0 from the [ Kudritzki et all 
d2000l) : ICiardullo et all d2002bl) : ICastro-Rodrfguez et all d2003h 
surveys, normalized to the same sampled volume. Thus, the ex- 
pected number of Ly-a galaxies are 65/15/6 in LPRX, 65/15/6 
in LPE, 66/16/6 in LPS and 86/20/8 in LSF field. 
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The blank field survey of ICastro-Ro dnguez e t alj d2003l) 
selected point-like emission line sources with same criteria 
as the current sample, and it provi des a lower inf e rred num- 
ber density of Ly-a th an either the ICiardullo et al.1 d2002bl) or 
KudritzkietaEI (2000) surveys. By statistically subtracting the 
contribution from the backg round population estimated from 
ICastro-Rodriguez et al.l d2003l) to the outer Virgo fields, we de- 
rive an upper limit to the number of ICPNe. The number of 
emission-line objects detected in the di fferent fields, and the ex- 
pected number of Ly-a objects from the ICastro-R odrigue z et alj 
(120031) survey brighter than m(5007)= 27.0 are 6/6, 14/6, 6/6, 
3/8 for LPRX, LPE, LPS, and LSF fields, respectively. The LPE 
field shows an excess of 8±4 emission-line objects over the ex- 
pected number of background galaxies in this field, while the 
LSF field shows a lower density than expected, by 5+2. We may 
therefore conclude that the emission-line objects detected in the 
LPRX, LPS, and LSF, fields are compatible with being all Ly-a 
background galaxies. The emission-line objects detected in LPE 
are compatible with containing a few ICPNe at the 2<x level, as 
in the SUB2/3 fields. 

An over-density of emission-line unresolved sources is de- 
tected at the RCN1 field position also; here we select 28 candi- 
dates using the color excess criteria in the CMD, using a deeper 
off-band image tha n in Paper I, and we expect a p opulation of 8 
Ly-a emitters from ICastro-Rodriguez e t al. (2003) scaled num- 
ber counts. 

On the basis of the previous results, the properties of the un- 
resolved emission line objects in these outer fields are those of 
the Ly-a population at z — 3.14. Within our sample, we find 
that the flux of the brightest emission sources may vary by a fac- 
tor 2 in differen t fields, i.e. from 26.0 in LSF field to 26.7 in 
SUB2 field. Th e Kudritzki et all (l2000l)lciar dullo et a i] J2002bl) . 
ICastro-Rodriguez et al.l d2003l) and lGronwall et alj ( 2007b sam- 
ples show similar field-to-field variation of the Ly-a LF bright 
end. Since these samples come from surveys that may adopt dif- 
ferent selection criteria, we believe that our current results on 
the field-to-field variations of the Ly-a LF in the six Virgo fields 
is more robust because it is based on homogeneously selected 
samples of point-like objects, with deep off-band control images. 
The variation of the Ly-a properties may be related to th e cos- 
mic variance: recent observations by Ouc hi et al] d2003l) show 
that Ly-a galaxies are already clustered at z=4.86, which implies 
even greater clustering at za 3. 

5. ICPN number density and ICL surface brightness 
in the Virgo cluster 

On the basis of the statistical subtraction of the contribution from 
Ly-a emitters at z = 3.14, we can conclude that: 

- In the SUB2 and SUB3 fields, we detected 4 ICPNe in both 
[OIII] and Ha. These are all brighter than m(5007)= 27.0, 
i.e. within 0.5 magnitudes of the PNLF bright cut-off at the 
distance of the Virgo cluster; 

- the number densities and LFs of the emission-line objects 
detected in the outer Virgo fields - LPRX, LPS, and LSF 
- are statistically consistent with those of the population of 
Ly-a background galaxies; 

- we detect over-densities in LPE (at the 2cr level) and in 
RCN1 (larger than 3cr) Virgo fields; 

On the basis of the number of emission line objects detected 
in the LPRX, LPS, LSF fields, we can set upper limits to the 
ICPNe population in each field. We statistically subtract the most 



conservative number of Ly-a emitt e rs exp ected in each field 
based on the ICastro-Rodriguez et al.l d2003l) survey (see discus- 
sion in Section |4~3| ). estimate the error-bars from Poisson statis- 
tics, and compute the lcr upper limit for the number of ICPNe 
independently for each field. For LPRX, LPS, and LSF, the ex- 
pected number of ICPNe brighter than m(5007)= 27.0 are 0+^, 
0^2 4, and -5+|*; for LPE and RCN1 the number of ICPNe is 
8+? g and 20+1 respectively. 



5.1. From ICPN number density to ICL surface brightness 

Since PNe closely fo llow the distribution of starlight in galaxies 
dCoccato et all I2009I) . the bolometric luminosity of the parent 
stars of the PN population can be obtained as L/, / = npn/a, 
where rippj is the number of PNs in each field and a is the 
luminosity-specific PN density. We follow the procedure de- 
tailed below: 

- the ripN for a given field is scaled by a factor: 



A = 



j^f +I PNLF(m)dm 

jm ltm p NLF ( nt)dm 



(1) 



where PNLF(m) is th e analytic expression for the PNLF 
dCiardullo et all Il989l) . M* and m* denote the absolute and 
apparent magnitude of its bright cutoff, respectively, and mu m 
is the m(5007) limiting magnitude in each field. This scaling 
ensures that we account for all PNs within 1 mag of M* . We 
then use the corresponding value a\ o to infer the amount of 
bolometric luminosit y in our fields. 

- iDohertv et"aT] d2009l) computed the a 2 .s for the halo of M87 
and the ICL directly. The value of a 2 . 5 ,MS7 is 3.1 x 10" 9 PN 
L~'; for the ICL a 2 . s.icL = 7.2 x 10" 9 ' PN L '; we refer to 
IDohertv et al.1 (120091) for further details on their computation 
and uncertainties. 

- From the analytical expression of the PNLF, we can scale the 
value of the luminosity-specific PN density at 2.5 mag fainter 
then the PNLF bright cut off to 1 mag as a\ q - 0.25 X a 2 .5- 

- a parameter for the FCJ field - The re s ults o f the 
spectroscopic follow -up by Arnaboldi et al. (2004) and 
IDohertv et all d2009i) show that the halo of M87 is extended 
out to about 35 arcmin from the galaxy center. Within that 
radius, the M87 halo contributes 70% of the light as sam- 
pled by the number of PNe detected at the position of the 
FCJ field. For this field position, we use a value of a\ o = 
1.1 x 10~ 9 PN L" 1 which is the average value of a\ .o,M87 
and ctisjcL, weighted by the fractions of PNs that are bound 
to M87 ( 70%) and those unbound (30%). Because the FCJ 
field is at a distance of about 14 arcmin from the center of 
M87, we can check the derived surface brightness value at 
this field posit ion with the V band surface brightness mea- 
surements of Korm endv et al.1 d2009l) . 

- a parameter for the SUBC field - we use a\ o = 1.3 X 10~ 9 
PN L ', which is the average of ai.o,M87 and ai.ojcL- Here 
we assume that M84 has the same ai.o value as M87 and 
the fractions of bound and unbound PNs in this photometric 
sample is 50% each. 

- a parameter for the outer fields - For the fields at larger dis- 
tances, we use ai.ojd - 4.0 x 10~ 9 PN L '. 

- Errors on the inferred luminosities and surface brightness are 
estimated from Poisson statistics. 

We compute the bolometric luminosities at each field position 
from the number of observed ICPNe, or their upper limit, and 
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Table 4. Npm within m * +1, B band luminosity density and sur- 
face brightness values in the Virgo fields. In LPC, LPRX, LPS 
and LSF we list the upper limits. 



Field 


N PN 


1 








10 6 L s o arcmirr 2 


mag arcsecr 2 


Core 


17.0 


3.1+0.7 


29.0+0.25 


FCJ 


15.0 


12.2+3.9 


27.5±0.25 


LPC 


< 1.0 


< 0.2 


> 32.0 


SUBC 


14.0 


4.0+0.7 


28.71±0.3 


LPRX 


< 2.4 


< 0.4 


> 31.1 


LPE 


8.0 


1.5±0.5 


29.8± 0.35 


LPS 


< 2.4 


< 0.4 


> 31.1 


LSF 


< 2.4 


< 0.4 


> 31.1 


RCN1 


13.0 


2.4±0.6 


29.3+ 0.3 


SUB2 


1.0 


0.2+0.2 


31.8+ 1.0 


SUB3 


3.0 


0.7+0.3 


30.6+ 0.6 



Table 5. Surface brightness measurements in the Feldmeie r et al.l 
(2004a) fields. 



Field 


UlCL.V 


fllCL.B 


B-V 




mag arcsec -2 


mag arcsec" 2 




F04-2 


27.4 


28.3 


0.96 


F04-3 


26.5 


27.5 


1.02 


F04-4 


27.4 


28.4 


1.02 


F04-5 


27.3 


28.3 


1.02 


F04-6 


28.1 


29.1 


0.96 


F04-7 


28.4 


29.4 


1.02 


F04-8 









sition, and 30% of the sample is made of free-floating stars, i.e. 
the true ICL component. The surface brightness of the ICL at the 
position of the FCJ field is then [i B = 28.8 mag arcsec -2 . 

5.3. How well do ICPNe follow light 

Mih os et alJ d2009) investigated the correlation between the 
ICPN distribution and the diffuse light in the Virgo core 
(Mih os et all l2005h . Their results indicate that there is a cor- 
respondence on large scale (~ 100 kpc), which is more ro- 
bust when the ICPN catalogs have a lower fraction of back- 
ground cojitaminants, as in t he case of the SUBC field (Paper 
II and Okamu ra et al. I (I2002ln . or in the higher surface bright- 
ness regions within the extended galaxy halos. A possible 
source of scatter is also a dependence of the luminosity- s pecific 
PN parameter on s tellar populations (Ciardu lfo et all 12005; 
iBuzzoniet a0.l2006b . 

The procedure adopted in this work maps the over-density of 
point-like emitters with respect to a statistically averaged back- 
ground population, in fields whose area is about 150 x 150 kpc 2 . 
As described in Section [5Tl the number of PNe from each over- 
density is mapped into a Bolometric luminosity according to 
the most appropriate luminosity-specific P N density parameter 
value for that field. Our comparison with the lMihos et ail (2009) 
results indicate that the luminosities computed in this work are 
robust with respect to the variety of effects that can act to wash 
out the correlation between PN number density distribution and 
light. 



the a parameter. Usi ng the bolometric correction to the B band 
dBuzzoni et al 1 120061) and (B-V) colors, we derive the B band 
luminosities and surface brightnesses. The adopted B-V colors 
are 0.9 for the fields at large distances, and (B-V) = 1.02 for 
those fields near either M87 and M84. In Table g] for each field 
we list the number of PNe within m* + 1 , Npn, the resulting mean 
values of the B-band luminosity, S, and surface brightness fig- 

5.2. ICL meas urements in the Virgo cl uster from previous 
work: the Feldm eier et al\ \2004a\) sample 

Se veral regions in the Virgo core and aro und M49 were surveyed 
by IFeldmeier etail dl998l 120031, l2004al) for ICPNe, covering a 
total area of 0.89 deg 2 . The positions of these fields are shown in 
Fig. [TJ and a summary of their fiy measurements is provided in 
Table|5] They find an excess of emission line sources in all fields, 
at about 5 times the density of line emitters in blank fields, but 
for Field F04-8, where the emission-line object number density 
is consistent with the galaxy background population. To compare 
these measurements with ours, we take their fi v and add a (B-V) 
color term determined from the color of the continuum light of 
the nearest bright galaxy; the computed values are shown in col. 
#4, in Table 

The field with the brightest surface brightness measurement 
is F04-3 (FCJ in this work; see also Paper I) at 60 kpc from 
the center of M87. This field is in fact within the M87 halo, as 
shown by the deep photometry of the Virgo core bylMihos et al.l 
d2005l) and the PN kinematics (lArnaboldi et all |2004|) . The n B 
measurement for F04-3 agrees with the independently derived 
fj.B for FCJ in Table [4] a nd with the s urface brightness profile of 
M87 measured by iKormendv et al.l d2009l) at the FCJ distance 
from the center of M87. The spectroscopic follow-up indicates 
that the M87 halo contributes 70% of the light at this field po- 



6. Discussion 

6.1. The diffuse light in the Virgo cluster 

Several fields outside the Virgo core region contain a population 
of emission sources consistent with background Ly-a objects. 
Clear over-density of emission line objects of at least a factor 
5 with respect to the averaged Ly-q emitters are detected in th e 
Virgo core region, around M49 (IFeldmeier et all I2003L l2004al) . 
and in RCN1 and LPE, at a lower <x level. 

The ICL and the extended halos from bright galaxies- The 
diffuse light observed in the core of a galaxy cluster contains 
several luminous stellar components that add up along the line- 
of-sight (LOS) to the cluster center: the extended faint halos of 
the brightest galaxies and the ICL contribution, defined as the 
light coming from stars not bound to individual galaxies. When 
computing the ICL fraction in the Virgo core, the surface bright- 
ness measurements must be corrected for the fraction of stars 
boun d to the extended halos of in dividual galaxies. The results 
from lArnaboldi etall d2004l) and iDohertv et all d2009l) on the 
LOS velocity distribution of the PNe indicate that the halos of 
M84, M86 and M87 are extended, out to 150 kpc in the case of 
M87. 

When we select only true ICPNe, we measure a surface 
brightness for the ICL of about fis = 28.8 mag arcsec -2 in FCJ 
and 29.0 - 29.5 in Core/SUBC. These surface brightness val- 
ues a re similar to those inferred from the detection of IC RGB 
stars dFereuson et all [19981; iDurrell et all 120021; Iwilliams et all 
12007b . 

The results from the spectroscopic follow-up imply that the 
measurements of the ICL based on a surface brightness thresh- 
old, e.g. jj v > 26.5, will lead to an overestimate of the ICL, since 
it includes the light contributions from extended galaxy halos. 
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6.2. On the elongated distribution of the ICL 

The iFeldmeier et al.l d200 4a)'s "smooth" elongated distribution 
for the ICL, of about 4 Mpc along the LOS to the Virgo core, 
was derived from the following observational results: 

- the relative brighter cut-off of the IC PNLF in sub-clump 
A, around M87, vs. sub-clump B, around M49, placing the 
sub-clump A at 4 Mpc in front of sub-clump B, 

- same number density of ICPNe in these fields. 

The evidence for a substantial depth along the LOS is weakened 
by the results that som e brightening of the PNLF in the M87 
halo may be intrinsic (Arna boldi et al.L 120081) . The presence of 
an extended halo a round M49, which i s to be expected given 
the results for M87 ( Doherty_etaLl 2009), may explain the same 
number density of PNe around M49 as in M87, because the M49 
halo reaches similar surface brightness as the M87 one at these 
fields' relative distances from the galaxy center. 

The upper limits in the LPC, LPS and L PRX fields from 
this w ork, and no-PN detection in F04-8 (Feld meier et all 
2003) show that there is not a uniform distribution of ICL all 
across the region fr om sub-clump A-M87 to subclump B-M49 
(Binggeli et allll987t) . The ICL drops below the detection limit 
in the region of the sky between these two Virgo sub-clumps. 
The current results indicate that the ICL is mostly associated 
with, and confined to, these high density regions. 

6.3. The ICL surface brightness radial profile 

In Fig.|6]we now plot the surface brightness values from Table|4] 
and [5] as functio n of the field distan ces from M87, for an easier 
comparison with lMihos et alJd2005l) and lKormendv et all 12009) 
photometry. 

Properties of the ICL distribution in the Virgo cluster - The 
ICPN survey in the Virgo cluster indicates that the ICL is not ho- 
mogeneously distributed in the core and at larger cluster radii. In 
the cor e, we have an upper limit for the LPC field, and in Field 
F04-8 (IFeldmeier et ali,l2003ll2004al) . Furthermore the deep im- 
age of the central region of the Vir go cluster shows the presence 
of complex structures at all scales (Mi hos et all [2005 ). Both re- 
sults indicate that the distribution of the diffuse light in the cen- 
tral region is inhomogeneous, further suppor t ing th e evidence 
that t he core is no t relaxe d, see lBinggeli et al.l (fl993). Paper III, 
and lDohertvetal.1 d2009l) . 

The comprehensive summary of all surface brightness mea- 
surements in Figure [6] indicates that most of the diffuse light 
is detected in fields located in the core, within a distance of 
about 80 arcmin from M87. Outside this region, the mean sur- 
face brightness decreases sharply, and the ICL is present in iso- 
lated pointings, e.g. RCN1 and LPE. 

Diffuse stellar light is also measured in "sub-structures", 
around M49 and i n the M60/M59 sub-gro up. The fields F04- 
2 and F04-6 from Feldmeie r et al.l d2004a) are situated in the 
outer regions of M49, at about 150 kpc from the galaxy center. 
These fields may contain PNe from the halo of M49 and the ICL 
component, w hich may have form ed within sub-clump B of the 
Virgo cluster dBinggeli et all Il987l) . Because the spectroscopic 
follow-up in not available yet for these PN candidates, we can- 
not quantify the fraction of light in the M49 halo and ICL for 
these fields. 

Using a technique similar to iDurrell et al.l d2002l) and 
IWilhams eta!1d2007l)jYan et alJd2008l) detected serendipitously 
an excess of point-like objects in two ACS fields at 7 and 15 ar- 
cmin from M60. The number density of the excess field popula- 



tion shows a decreasing gradient with increasing distance from 
M60. This population starts at m/775 = 26.2 and the average sur- 
face brightness is 26.6 mag arcsec 2 in the I band in the field 
at 7 arcmin from M60. These va lues are very similar to those 
measured by IDurrell et all d20 02) observations of the IC RGB 
in the Virgo core. lYan et all d2008l) argue that they detected a 
stellar population associated with the M60/M59 pair. This result 
may indicate that the M60/M59 sub-group is yet another region 
of the Virgo cluster where the diffuse halos and Intra-group-light 
are being assembled. 

The current result on slope of the radial profile of the ICL 
is somewhat differen t from the slow decline in ICPN number 
density described by IFeldmeier et all d2004al) . A steeper radial 
ICPN density profile is now derived because the current survey 
probes areas between the sub-clump A and sub-clump B in the 
Virgo cluster, where only upper limits are set. 

6.4. Fraction of the ICL light vs. galaxy light in the Virgo core 

When we average the measurements of the ICL in the core re- 
gion defined as the circular area of 80 arcmin radiu^] centred 
on the Virgo cluster center given by Binggeli et al. (1987), 
we obtain a measurement of fig = 29.3 mag arcsec~ 2 . The 
ICL luminosity in t his area is L q jcl = 39.2 x 10 9 L O; g. From 
iBinggeli et al.l dl987l) the cumulative luminosity from all galaxy 
types within 1.3 deg = 80 arcmin radius of the Virgo cluster cen- 
ter is L^gai = 58.16 x 10 10 L Qi b. Our best estimate of the fraction 
of ICL, defined as the un-bound stellar component in the Virgo 
cluster core, amounts to ~ 7% of the total light in Virgo cluster 
galaxies. 

T he current fraction is lo wer than the recent measurement 
from Feldmei er et alj d2004al) . using a very similar technique. 
This is not unexpected: we have corrected the surface brightness 
measurements in the inner Virgo core fields by subtracting a sig- 
nificant fraction of PN bound to bright galaxies halos, which is 
between 70% and 50%, of the total Npn measured in those fields. 
We think that the current (conservative) ICL fraction is a robust 
estimate of the light from true-unbound stars within 80 arcmin 
centred on M87. 

6.5. Implications for the formation of the ICL in the Virgo 
cluster 

We can summarize the current results as follows: 

- The ICL in Virgo is concentrated in the core, with Hicl,b in 
the range 28.8 - 30 mag arcsec~ 2 . 

- In several fields at R > 80' = 340 kpc we can set only upper 
limits to the ICL surface brightness. 

- At R > 80', the ICL signal is confined to isolated fields 
(RCN1,LPE) 

- Diffuse stellar light defined as extended halo plus ICL is 
present in other Virgo sub-structures (central region of sub- 
clump B, M60/M59 sub-group). 

The distribution of the ICL in the core and at larger scales in 
the Virgo cluster shows that the ICL is more centrally concen- 
trated than the galaxi es, see Fig. [6] This r esult is consistent with 
the measurements of I Zibettie t al. (2005) and the predictions of 
iMurante et ail d2004|]2007l) . 



4 Because this circular area includes regions with only upper lim- 
its, our ave r age s urface brightness is lower than that measured by 
iMihos e t al. ( 200^) in the region between M87 and M86. 
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100 150 200 
R (arcmin) 

Fig. 6. Surface brightness measurement of diffuse light in the 
Virgo fields (points) compared with the surface brightness pro- 
file of the Virgo galaxies averaged in annuli (lines); radial dis- 
tances are computed with respect to M87. The green line rep- 
resents the ra dial surface brightnes s profile from light in Virgo 
galaxies from Bing geli et al.ldl987l) . The dotted-dashed and dou- 
ble dotted-deshed lines correspond to the surface brightness 
profile associated with giants and dwarf galaxies, respectively. 
The full blue dots show the surface brightness measurements 
in the FCJ, Core, SUBC and LPC fields. The open circles at 
the FCJ, SUBC field position indicates the ICL surface bright- 
ness computed from the ICPNe not bound to galaxy halos. The 
triangles represent t he surface brightness of the ICL based on 
IC RGB star counts dFerguson et al.L I1998I; iDurrell et all l2002t 
I Williams et al.Ll2007h . The full red dots show the surface bright- 
ness measurements in the RCN1, LPE, SUB2, and SUB3 fields. 
The arrows indicate the upper limits for the measurements in 
LPRX, LPS and LSF The magenta open squares indicate the 
surface brightness average values fig at 15, 50 and 240 ar- 
cmin computed from th e measurements listed in Table 5, from 
iFeldmeier et alj d2004al) data; the measurements at 240 arcmin 
(F04-2 and F04-6) are close to M49, see Figure [Q Radial dis- 
tances are computed with respect to M87. 



The ICL distribution shows a high degree of substructures 
both in the Virgo core and at larger radii. This result is sup- 
porte d by deep surface brightness photometry ( Mihos ^t al.L 
120051) . and ICPN number counts. The presence of sub-structures 
is predi cted by cosmological s imulations of structu r e for- 
mation (Napolitan o et all l2003t ISommer- Larsen et all [2005; 
iMurante et all 120071) . The detection of extended halos and 
intracluster/intra-group light in Virgo sub-structures supports a 
hierarchical formati on mechanism for the ICL as described by 
iRudick et all (2006): we are observing the ICL component in 
sub-structures, which may then become completely unbound as 
these sub-structures fall into the densest part of the cluster. 

Current data are consistent with most ICL being associated 
to the formation history of the brightest cluster galaxies, which 
explains the central concentration of the ICL in the core and to 
sub-structures. Tidal effects are indeed at work at larger dis- 
tances from the center and they can generate some ICL. The 
ICPNe detected in RCN1, LPE and SUB2/3 would likely come 
from such processes, but the amount of light detected in these 
fields is small compared with that sampled close to the big ellip- 
ticals, as supported by the observational results of this work. 



7. The Local Mean Luminous density 

From the lack of PN detections in several of our fields in the 
Virgo cluster, we can derive a limit on the luminosity density in 
an as sumed homogeneous distribution of stars in the local vol- 
ume. Fukugita et alJd 1998b compiled measurements of the mean 
luminous density contributed by high surface brightness galaxies 
at moderately low redshift, and adopted eo = 2x IQ & hLg^M pc~ 3 . 
In what follow, we make the assumption that the mean luminous 
density is distributed uniformly in the local volume and compute 
the number of PNe associated with this diffuse component in the 
surveyed area of the Virgo outer fields and at the [OIII] limiting 
magnitude m/,„, (5007) = 27.0 

Each volume element along the LOS to the Virgo fields con- 
tribute the following number of PNe 



dN PN (r) = a(r) x dUf) 



(2) 



where L(r) is the luminosity associated with the volumetric lu- 
minous density 



dL(r) = eo x S x dr 



eoQ, 



206265 2 



r 2 dr 



(3) 



e = 2 x 10*hL BtO Mpc~ 3 (see iFukugita etafl Il998), Q. is the 
area in arcsec 2 , a(r) is the luminosity-specific planetary nebula 
number and predicts the number of PNe associated with the to- 
tal luminosity L(r) at radial distance r. Its value depends on the 
relative depth of the on-band survey (/% m ) with respects to the 
faint limit of the PNLF, e.g. the apparent magnitude of M* + 8 at 
the distance r. 

The total number of PNe associated with the volumetric den- 
sity is then given by the integral of eq. [2] 



N, 



PN,tol 



f N PN {r)dr= f an 

Jo Jo 

I a 

Jo 



(r)L T (r)dr 



(4) 



206265 2 



(rXdr 



and r ma x is set by the survey depth, and it is the distance at which 
the apparent magnitude of the PNLF bright cut-off m* is equal 
tom, im . 

The luminosity-specific planetary nebula number a - The 
PN population (N/w) associated to a stellar population with lu- 
minosity Lj is : 



Npn = 



M*+8 



PNLF(M)dM = 



(5) 



ke 
Jm" 



,0.307 M 



(1.0 -e 



(3*(M*-M)) 



)dM = aL T 



where M* = -4.5 is the absolute magnitude of the bright cut-off 
of the P NLF, and a is the lum inosity-specific PN number, a is 
given by iBuzzoni et alj d2006l) for stellar populations with dif- 
ferent ages and metallicities. The maximum theoretical value of 
a which gives the largest PN population for a given luminosity 
is a max = IPN x (1.85 x 1O 6 L ) _1 . From eq.flwe see that the 
PNLF extends over 8 magnitudes, but one sees smaller magni- 
tude ranges as distance increases. The a parameter must then be 
corrected using the double-exponential PNLF analytical formula 
in eq.|5]for the appropriate faint cut-off magnitude set by the the 
minimum between m;, m and the apparent magnitude correspond- 
ing to M* + 8 at distance r. 

Predicted number of PNe - We use eq. [4] to compute the 
number of PNe we would predict if the Fukugita et al. ( 1998) 
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density was homogeneously distributed. The number of PNe that 
we would detect by surveying 1.75 deg 2 (total area covered by 
the L* fields plus SUB2 and SUB3 fields) down to m Vm = 27.0, 
is NpNjot = 4.6 with an expected m(5007) magnitude distribu- 
tion that peaks at ~ 26.0. We conclude that we cannot set limits 
on the homogeneously distributed luminous density in the local 
universe, and surveys of several hundreds square degrees are re- 
quired to detect a statistical significant numbers of the brighter 
PNe. 



8. Conclusions 

We have surveyed several regions of the Virgo cluster located at 
different distances from M87, covering a total area of 3.3 square 
degrees. These outer fields in the Virgo cluster were imaged 
through off- and on-band filters designed for detecting PNe at the 
distance of the Virgo cluster. The result of this survey provides 
an homogeneous large scale study of the diffuse stellar light in 
Virgo. 

The results from the current study indicate that the diffuse 
stellar component is present in the central region of the cluster, 
within a radial distance of 80 arcmin from the cluster center. 
The intracluster light in the core has a surface brightness of /ig = 
28.8 - 30 mag arcsec -2 , and on average it amounts to ~ 7% 
of the total galaxy light in this region. The distribution of the 
ICL in the cluster central region is not homogeneous, with areas 
where only upper limits on the number of PNe can be set. This 
non-uniformity indicates that the ICL is not relaxed in the core. 

At a distances larger than 80 arcmin from M87, the surface 
brightness profile of the diffuse light has a sharp cut-off of about 
2 magnitudes. At these larger radii the detection of ICL is con- 
fined in single fields, e.g. RCN1, LPE and SUB 2/3, or associated 
with local substructures, e.g. sub-clump B, the M60/M59 sub- 
group. The observations indicate that the ICL in Virgo is more 
centrally concentrated than cluster galaxies, and does not sup- 
port an homogeneous elongat ed distribution of the ICL similar to 
the la te-type Virgo galaxies jYasuda et all 119971: ISolanes et al.L 
l2002h . 

The concentration of the ICL toward the cluster high density 
regions links the formation of the ICL with the formation and 
evolution of the most lum inous cluster galaxies dRudick et al.L 
2006: lMurante et aUl2007l) . 
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Appendix A: Catalogs of the point-like emission 
line objects detected in the current survey 

Table A.l. Emission-line objects detected in the LPRX field 



Name a(72000) (5(72000) m(5007) 



RCN1-1 


12: 


:25: 


49 


.69 


14: 


:22: 


:08 


.86 


25 


,478 


RCN1-2 


12 


:27: 


11 


.53 


14: 


:06 


:36 


.37 


25 


,913 


RCN1-3 


12: 


:27: 


:04 


.63 


14: 


:14: 


:55 


.15 


25 


,915 


RCN1-4 


12: 


:26: 


01 


.43 


14: 


19 


:13 


.80 


25 


.926 


RCN1-5 


12: 


:26: 


11 


.00 


14: 


:08 


:49 


.73 


25 


.926 


RCN1-6 


12: 


:26: 


:30 


.36 


14: 


:07 


:31 


.63 


26 


.194 


RCN1-7 


12 


:26: 


:40 


.40 


14: 


:08 


:35 


.24 


26 


.240 


RCN1-8 


12 


:27: 


:12 


.32 


14: 


:01 


:27 


.14 


26 


.249 


RCN1-9 


12 


:26: 


:30 


.01 


14: 


:17: 


:09 


.80 


26 


,295 


RCN1-10 


12 


:25: 


:34 


.38 


13: 


:58 


:09 


.08 


26 


,297 


RCN1-11 


12. 


26: 


:05 


.27 


14: 


:11 


:14 


.95 


26 


,311 


RCN1-12 


12 


25: 


:35 


.01 


13: 


:59 


:42 


.98 


26 


,312 


RCN1-13 


12 


21: 


:13 


.55 


14: 


:17 


:48 


.79 


26 


.343 


RCN1-14 


12 


26: 


59 


.59 


14: 


:09: 


:56 


.40 


26 


,348 


RCN1-15 


12 


21: 


08 


.51 


14: 


:08 


:29 


.50 


26 


,355 


RCN1-16 


12: 


26: 


59 


.25 


13: 


:53 


:42 


.86 


26 


.359 


RCN1-17 


12: 


26: 


28 


.77 


14: 


:08 


:34 


.98 


26 


,394 


RCN1-18 


12 


26: 


33 


.17 


14: 


:16 


:20 


.35 


26 


.489 


RCN1-19 


12 


26: 


:28 


.15 


14: 


:08: 


:04 


.62 


26 


,490 


RCN1-20 


12 


25: 


:55 


.82 


14: 


:09 


:08 


.78 


26 


,572 


RCN1-21 


12 


21: 


:04 


.90 


14: 


:10: 


:43 


.36 


26 


.602 


RCN1-22 


12: 


25: 


:54 


.25 


13: 


:55: 


:17 


.77 


26 


.657 


RCN1-23 


12 


21: 


15 


.57 


14: 


:23 


:12 


.65 


26 


.674 


RCN1-24 


12 


21: 


:07 


.64 


14: 


.16 


:10 


.91 


26 


,681 


RCN1-25 


12: 


25: 


:27 


.23 


14: 


:11 


:33 


.21 


26 


.682 


RCN1-26 


12: 


21: 


:07 


.01 


13: 


:59 


:56 


.88 


26 


.686 


RCN1-27 


12: 


26: 


:34 


.64 


14: 


:11 


:55 


.75 


26 


,689 


RCN1-28 


12: 


26: 


:09 


.25 


13: 


:57: 


:22 


.98 


26 


,688 



Table A.2. Emission-line objects detected in the LPRX field 



Name 


a(72000) 


5(72000) 


m(5007) 


LPRX-1 


12:27:20.377 


9:17:37.14 


26.03 


LPRX-2 


12:26:32.996 


9:23:48.76 


26.42 


LPRX-3 


12:28:11.142 


9:25:14.66 


26.47 


LPRX-4 


12:27:25.608 


9:16:14.40 


26.53 


LPRX-5 


12:27:48.012 


9:14:26.81 


26.78 


LPRX-6 


12:26:42.855 


9:16:12.58 


26.94 


LPRX-7 


12:27:25.148 


9:12:53.31 


27.19 


LPRX-8 


12:27:31.843 


9:37:10.15 


27.33 


LPRX-9 


12:26:40.641 


9:28:57.08 


27.36 



Table A.4. Emission-line objects in the LPE field 



Name a(72000) 5(72000) m(5007) 



LPE-1 


12: 


17 


:58 


,471 


13 


:40 


:07 


.24 


26.08 


LPE-2 


12: 


17 


:32 


.446 


13 


:48 


:31 


.87 


26.15 


LPE-3 


12: 


17 


:39 


,261 


13 


:32 


:27 


.24 


26.19 


LPE-4 


12: 


18 


:04 


.059 


13 


:43 


:38 


.15 


26.19 


LPE-5 


12: 


17: 


:59 


.871 


13 


:43 


:18 


.95 


26.40 


LPE-6 


12: 


17 


:17 


,932 


13 


:28 


:25 


66 


26.43 


LPE-7 


12: 


16 


:04 


.795 


13 


:40 


:59 


.73 


26.44 


LPE-8 


12: 


16 


:36 


.738 


13 


:45 


:45 


.07 


26.52 


LPE-9 


12: 


16 


:02 


.328 


13 


:31 


:03 


.16 


26.71 


LPE- 10 


12: 


17 


:01 


.922 


13 


:56 


:48 


.82 


26.8 


LPE- 11 


12: 


18 


:04 


.406 


13 


:43 


:05 


.79 


26.8 


LPE- 12 


12: 


16 


:54 


,828 


13 


:58 


22. 


.60 


26.9 


LPE- 13 


12: 


18 


:14 


.529 


13 


:39 


:13 


.96 


27.0 


LPE- 14 


12: 


18 


:07 


.388 


13 


:52 


:02 


.52 


27.0 


LPE- 15 


12: 


16 


:53 


,685 


13 


:29 


:48 


.63 


27.1 


LPE- 16 


12: 


16 


:48 


,095 


13 


:52 


:14 


.56 


27.2 


LPE- 17 


12: 


16 


:59 


.375 


13 


:44 


:38 


.69 


27.3 


LPE- 18 


12: 


16 


:41 


.977 


13 


:34 


:25 


,83 


27.3 


LPE- 19 


12: 


18 


:09 


.217 


13 


:50 


:04 


,31 


27.4 


LPE-20 


12: 


18 


:11. 


.198 


13 


:57 


:59 


.73 


27.6 



Table A.5. Emission-line objects detected in the LSF field 



Name a(72000) 5(72000) m(5007) 



LSF-1 


12: 


:37: 


:42 


.215 


12: 


15 


■A3. 


.10 


25 


99 


LSF-2 


12: 


:39: 


:13 


.235 


12: 


:14 


:38. 


.27 


26. 


.83 


LSF-3 


12: 


:38 


:07 


.374 


12: 


:07 


:34. 


,85 


27 


.08 


LSF-4 


12: 


:39: 


:13 


,669 


12: 


01, 


:46. 


,78 


27 


.17 


LSF-5 


12: 


:37: 


:42 


,473 


12: 


:11 


:20. 


.97 


27. 


.35 


LSF-6 


12: 


:39 


:20 


.662 


12: 


16 


:37. 


.00 


27. 


.36 


LSF-7 


12: 


:39 


:09 


.730 


12: 


10 


:50 


.76 


27 


.37 


LSF-8 


12: 


:39 


:13 


.045 


12: 


:11, 


:06. 


,79 


27 


.66 


LSF-9 


12: 


:37 


:46 


.282 


11: 


:55, 


:14. 


,22 


27 


.70 


LSF- 10 


12: 


:38 


:47 


.034 


12: 


18 


:00. 


90 


27 


.70 


LSF- 11 


12: 


:37 


:28 


,577 


12: 


:17: 


:48 


,96 


27 


.80 


LSF- 12 


12: 


:37 


:46 


,895 


12: 


15 


:24 


,38 


27 


.91 


LSF-13 


12: 


:37: 


:38 


,274 


12: 


:20: 


:31. 


,31 


27 


.93 


LSF- 14 


12: 


:39 


:15 


.965 


11: 


:55: 


:43. 


,51 


27 


.93 


LSF- 15 


12: 


:37 


:30 


.846 


12: 


19 


:05 


18 


27 


.97 


LSF- 16 


12: 


:39: 


:27 


.344 


12: 


:20 


:25. 


,23 


27, 


.99 


LSF- 17 


12: 


:39: 


:20 


.736 


12: 


:05, 


:44, 


.44 


28 


.02 


LSF- 18 


12: 


:37, 


:29 


.078 


12: 


:22: 


:37. 


17 


28 


.03 


LSF- 19 


12: 


:38 


:05 


.800 


12: 


01 


:34 


.96 


28 


.03 


LSF-20 


12: 


:38 


:41 


,080 


11: 


:57, 


:07 


,98 


28 


.04 


LSF-2 1 


12: 


:38 


:00 


,317 


12: 


21: 


:35. 


,34 


28 


.06 


LSF-22 


12: 


:38: 


:37 


.500 


12: 


:05 


:56. 


12 


28, 


09 


LSF-23 


12: 


:39 


:23 


.081 


12: 


:12 


:33, 


,49 


28 


.15 



Table A.3. Emission-line objects detected in the LPS field 



Name 


a(72000) 


5(72000) 


m(5007) 


LPS-1 


12:27:25.204 


10:51:38.34 


26.21 


LPS-2 


12:26:38.216 


10:35:49.34 


26.36 


LPS-3 


12:27:34.998 


10:26:46.14 


26.38 


LPS-4 


12:27:10.849 


10:43:44.31 


26.75 


LPS-5 


12:27:33.526 


10:44:30.01 


26.98 


LPS-6 


12:25:59.211 


10:44:33.17 


27.01 


LPS-7 


12:27:38.171 


10:53:50.34 


27.10 



N. Castro-Rodriguez et al.: Intracluster Light in the Virgo Cluster 
Table A.6. Emission-line objects detected in the SUB2 field. Objects with (*) are ICPNe 



Name 


a(72000) 


(5(72000) 


m(5007) 






SUB 2-1 


12:23:22.524 


14:03:21.68 


26.30 


— 


— 


SUB2-2 


12:24:23.113 


14:20:50.64 


27.01 


— 


30.22 


SUB2-3 


12:24:45.272 


14:18:52.58 


27.18 


— 


29.91 


SUB2-4 


12:24:35.834 


14:05:36.93 


27.19 


— 


29.26 


SUB2-5 


12:24:45.103 


14:16:12.88 


27.34 


— 


27.56 


SUB2-6 


12:24:46.116 


14:13:04.52 


27.46 






SUB2-7* 


12:24:15.861 


14:04:17.51 


27.54 


25.81 


28.33 


SUB2-8 


12:25:01.022 


14:02:10.35 


27.58 




28.53 


SUB2-9 


12:23:39.559 


14:18:15.09 


27.59 




31.73 


SUB2-10 


12:24:35.553 


13:55:21.40 


27.75 






SUB2-11 


12:23:59.558 


14:13:02.89 


27.85 






SUB2-12 


12:23:04.010 


14:06:12.86 


27.88 






SUB2-13 


12:23:25.975 


14:05:48.30 


27.90 







Table A.7. Emission-line objects detected in the SUB3 field. Objects with (*) are ICPNe 



Name g(72000) (5(72000) m(5007) m Ha m (v+R) 



SUB3-1 


12: 


:24 


:21.128 


15: 


:32. 


: 10.67 


24.57 






SUB3-2 


12: 


:25 


:02.814 


15: 


:16: 


:42.17 


25.75 






SUB3-3 


12: 


:25 


:05.471 


15: 


:16 


:30.21 


26.13 






SUB3-4 


12: 


:25 


:05.812 


15: 


:16 


:28.84 


26.25 




28.10 


SUB3-5* 


12: 


:24 


:26.934 


15: 


:15 


:43.10 


26.31 


25.46 


27.30 


SUB3-6* 


12: 


:25. 


:05.369 


15: 


:23. 


: 04.51 


26.46 


25.46 


27.66 


SUB3-7* 


12: 


:23 


:32.742 


15: 


:16 


:58.09 


26.49 




27.06 


SUB3-8 


12: 


:25: 


:15.651 


15: 


:27: 


:13.45 


26.57 




26.86 


SUB3-9 


12: 


:25 


: 00.771 


15: 


:15: 


:52.94 


26.99 




26.72 


SUB3-10 


12: 


:24. 


:35.389 


15: 


:08 


:58.35 


27.52 




28.72 


SUB3-11 


12: 


:23 


:40.561 


15: 


:09: 


:28.04 


27.53 







